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Tracking Pollen for Asthma Alerts in Public Health DSS (Luvall) 



Earth Observations 



MODIS & NPOESS-VIIRS 
MOD9, MOD09GA, 
MOD09SPA, NDVI/ 
EVLSPA, MODLSTSPA, 

ASTER 

Burkard pollen samplers 
ASD-FR 

spectroradiometer 
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transport 


Asthma/Ml health data 
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Optimize hospital staffing 
Allergists alerted to pollen timing 


OGC Web Map Services 
Apache Web Sever 
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Value and 
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and society 


alth alerts 


Reduced medical costs 
Fewer inpatients 
More accurate diagnoses 


Inputs Outputs 


NASA and UNM, UA Partners 


Outcomes 


Impacts 


MDOH and CDC EPHT Systems; SYRIS 
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Pollen and Respiratory Disease: 
Wha t litt le is known£ • • 


Increase in mortality of 
these disorders: 


Cardiovascular disease 

CtlfHlIc Obstructive pulrnuriury disease 

Pneumonia 

Total 


flovneeae pollen concentrations [grains per rtf air) j 

<22 

22-77 

78-115 

>135 

Relative risk 

Relative risk (05% Oil 

Relative risk (95% 01) 

Relative risk (95% 01) 

1000 

1-015 11 -002-1- 029) 

1-012 (0-394-1-029) 

1-061 (1-038-1-0841 

1-000 

1-035 '1-053-1 -139) 

1-124 (1-069-1- 131) 

1-150 (1-079-1-2251 

1-000 

1104 |1- 049-1- 163) 

1-093 (1-023-1- 168) 

1-168 (1-077-1- 2661 

1-000 

1-019 fl-OlO-1-028) 

1-019 (1-00 8-1-0311 

1-043 (1-028-1-0581 


• High concentrations of pollen allergens have also been shown to occur in thoracic particles (<10 
microns in diameter) and respirable particles (<2-5 microns and these correlated well in time with 
airborne pollen concentrations. ... airborne pollen results in exposure of the lower airways and 
lung to pollen allergens. 

• The association between air pollution and the number of daily deaths may be related to the 
inflammatory potential of very small particles 

• ...suggests that high airborne pollen concentrations, which nowadays are mainly seen as triggers 
of allergic symptoms, may have far more serious effects than previously thought. 


2 Bert Brunekreef, Gerard Hoek, Paul Fischer, Frits Th M Spieksma. Relation between airborne pollen concentrations and 
daily cardiovascular and respiratory-disease mortality. Lancet Vol 355 (2000): 1517-8. 
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Juniper Species and Pollination Season 


Juniper Ashei (J a ) is mostly found to be distributed over Texas and 
Oklahoma and pollinates during December to January. Thus, the 
dispersion of juniper pollens during December-January is mostly 
restricted to J a type.. 


Juiper Pinchotti (J p ) is mostly distributed over Texas and pollinated 
during October-November. Thus, the dispersion of juniper pollens 
during this period is mostly restricted to J p type. 


Juniper monosperma (J m ) and Juniper scopulurum (J s ) are 

prevalent in New Mexico and pollinates during March-April period. 
Thus, the dispersion of juniper pollens during this period is mostly 
restricted to J m and J s type. 



transport 


■ 27 Jan 99, Jim Anderson in 
London, Ontario reported 
atmospheric Juniperus pollen - 
58 pollen grains/m 3 

■ T rajectories show that the 
source of this pollen was Texas 
population of Juniperus ashei 

■ Our Jan 26 forecast indicated 
that the “pollen has the potential 
to travel very long distances.” 
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NOAA AIR RESOURCES LABORATORY 


Trajectory Forecasts 

Farrmd Ti^cdotj Staving- 17 UTC 26 JAN 99 
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Pollen per cone 


Sonora 
Dallas 
San Marcos 
Junction 


472.000 pollen grains/cone 

402.000 pollen grains/cone 

374.000 pollen grains/cone 

363.000 pollen grains/cone 


269,946 to 946,646 
cones per tree 



1 60 to 1 500 
Trees/Ha 


^9 





Pollen Release 








itations of Pollen Sarnpljng . 



Lack of stations 

Count frequency & reporting lag time 

Different sampling instruments Rotorod Sampler/Burkard 
Spore Trap 

Only indentifiable pollen “grains” 

Expertise in counting/indentification 

Refusal to release sampling information-” We do not 
reveal the sources for our data for privacy and 
proprietary, competitive reasons. Some pollen counts are 
conducted privately, and are not meant to be broadcast 
to the public " 




n Timing 



Growing Degree Days - the average of the daily 
maximum and minimum temperatures compared to a 
base temperature, T base , (usually 10 ° C) 

Response to length of day 

Species differences 

Climate - Variability in Precipatation 

Weather 




Southern Rocky Mountain Juniper Woodland and Savanna 
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►ectral characteristics of male juniper canopies 

* ■ . ■ 

t-dtfrerent bud density levels ' 


Density Bud density 
level (g/m A 2) 


204.2 

190.0 

176.9 

164.9 

151.1 

136.2 
115.8 

92.9 

45.9 
0.0 
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remporal profiles of residual MODIS reflectances at the four study sites 


Many challenges 

Residual signals and reference baselines 

Landscape vs species level phenology & 
signals (disaggregate woody from herbaceous) 

Surface heterogeneity and spatial 
characterization of landscape 

Future sensors & fusion (Lidar, VIIRS, 
HyspIRI) 

BRDF & surface aerodynamics 
Modeling (vegetation dynamics, phenology) 




Summary Statistics for 2010 and 2011 Juniperus 
pinchotii pollen seasons at source 


Location 

Average daily 
concentration 
Pollen 
grains/m3 

Peak daily 
concentration 
Pollen 
grains/m3 

Date of 
peak 

Peak hourly 
concentration 
Pollen grains/m3 

Time of 
peak hour 

Date of 
peak hour 

2010 

Erick, OK 

337 

5,563 

25-Oct 

15,898 

10:00 AM 

25-Oct 

Sonora, TX 

286 

3,019 

25-Oct 

12,152 

10:00 AM 

25-Oct 

San Angelo, TX 

653 

5,542 

1-Nov 

10,195 

Noon 

1-Nov 

2011 

Erick, OK 

12 

214 

16-Oct 

800 

10:00 AM 

16-Oct 

Sonora, TX 

64 

428 

26-Oct 

2,422 

10:00 AM 

31 -Oct 

San Angelo, TX 

60 

493 

26-Oct 

2,747 

4:00 AM 

26-Oct 

Quanah, TX 

22 

190 

19-Oct 

703 

8:00 PM 

19-Oct 


Influence of preseason meteorological 
variables for Juniperus ashei 

► Start date 

► Significantly correlated with mean monthly temperature 
in December (r = 0.467, p = 0.038) and November 
rainfall (r = 0.468, p = 0.038) 

► Cumulative Season Total (CST) 

► Significantly correlated with mean maximum temperature 
in December (r = 0.4740, p = 0.035) 


Trajectory Wind Directions 



Southerly 

Northerly 

Within TX 

Other 

2007 to 2008 Pollen Year 


17% 

54% 

24% 

5% 

2008 to 2009 Pollen Year 


12% 

68% 

7% 

12% 

2009 to 2010 Pollen Year 


31% 

50% 

7% 

11% 

2010 to 2011 Pollen Year 


26% 

41% 

30% 

4% 

2011 to 2012 Pollen Year 


18% 


54% 


16% 


13% 
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Trajectory analysis indicates the pollen originated in 
southwest Texas approximately 8 hours earlier 




NOAA HYSPLIT MODEL 

Forward trajectories starting at 1600 UTC 25 Oct 1 0 

GDAS Meteorological Data 


This is not a NOAA product It was produced by a web user. 

Job I D: 349594 Job Start: Thu Feb 1 6 20:06:47 UTC 201 2 

Source 1 lat.: 30.26 Ion.: -101 .56 height: 500 m AG L 
T raj ectory D i recti on : Fo rward D u ratio n : 24 h rs 

Vertical Motion Calculation Method: Model Vertical Velocity 

Meteorology: 000 QZ 22 Oct 2010 - GDASl 


500 


Juniperus pinchotii concentrations at Sonora 
and San Angelo confirm the trajectory model 




Hygroscopic Weight Gain of Juniperus 

Pollen 


• Juniperus ashei, J. monosperma, 
and J. pinchotii were dusted on 
greased slides and exposed to a 
range of temperatures and 
humidity levels and weighed at 2 
hrs and 6 hrs. 

•V. ashei was weighed at 
temperatures 4° C, 15° C, and 
20° C at relative humidity levels; 
20%, 40% 50% 75% 86% 97%. 


•Weight was not significantly 
affected by temperature or time. 



Estimated weight per pollen grain (J. ashei) after 2 hrs across the 
range of relative humidity levels at 20° C 





ne Dust Simulations and Forecasts 

ity of Arizona • ’N, v ■ * . . 

nee & University of New Mexico 
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Phoenix dust storm - 7 June 2006 
Photo by Robb Schumacher Arizona Repu 


Department of Atmospheric 
Sciences 


http://www.atmo.arizona.edu/faculty/research/dust/dust.html 
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UA WRF 10-m wind forecast 


S. Nickovic et al., A model for prediction of desert dust cycle 
in the atmosphere, JGR 106, 18113-18129 (2001) . 

Yin et al., Modeling wind-blown desert dust in the 
southwestern United States for public health warning: A case 
study, Atmos. Environ. 39, 6243-6254 (2005). 

Yin et al.,The impact of using different land cover data on 
wind-blown desert dust modeling results in the southwestern 
United States Atmos. Environ., 41, 2214-2224 (2007). 


Adapted from Betterton ppt 
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Dust REgional Atmospheric 
Modeling (DREAM) 
system 

■ MM5 

■ WRF 



Pollen REgional 
Atmospheric Modeling 
(PREAM) system 



Model horizontal domain : Southwest US 
Model resolution : ~19 to 40 km 
Bins: 4, sized by pollen gain size distributions 
Boundary conditions: driven by the NCEP/WRF 
1 degree global forecasts used to refresh; initial 
conditions every 24 hours, boundary conditions 
every 6 hours 

Model Output : Every 3 hr maps of pollen surface 
concentrations (grains/m 3 ) out to 48 hrs. 


UA WRF 10-m wind forecast 

S. Nickovic et al., A model for prediction of desert dust cycle in the atmosphere, JGR 
106, 18113-18129 (2001) . 

Yin et al., Modeling wind-blown desert dust in the southwestern United States for 
public health warning: A case study, Atmos. Environ. 39, 6243-6254 (2005). 

Yin et al., The impact of using different land cover data on wind-blown desert dust 
modeling results in the southwestern United States Atmos. Environ., 41, 2214- 
2224 (2007). 








Applications Have Included: 

■ Dust Storms & Airborne Mineral Dust 
Concentrations in the Middle East, Africa and the 
Southwest US 

■ Pollen in Colorado, New Mexico & Texas 
■Volcanic Ash in the Mediterranean 

■ Soybean Rust in South America 

A new test: Forest fire ash and smoke plumes 
A proposed test: mold spores 


VI 



M 4-8 particle 

— • / . 


Horizontal 

advection 


Model predictions (72-h): 

Horizontal distribution 

• Surface concentration 

• Total column mass (dust loa 

• Wet, dry, total deposition 

• Meteorological variables | 

Vertical distribution 

• Concentration 

t m Mineral 

• Cross sections 

• Fixed point/time profiles 

Fixed point (selected sites/cities) 



Absorption 

Reflection 





■ Select Pollen of Interest 

■ Map Pollen Source 

■ Estimate Emission on Test Date 

■ Prepare Model 

• Insert Terrain & Pollen Aerodynamic Characteristics 

• Insert Source Emission 

• Insert Meteorology 

■ Simulate Downwind Pollen Dispersal 

■ Evaluate 


Pollen release potential Source Map/Mask 


(PRPSM_of_Jj) of a Juniper species “/” is calculated as: 

PRPSM_of_Ji = Tj x Mj x Hj x Cj x P; 


Where 

Tj = Number of Jj trees 
Mj = Male/Female ratio of Jj 
H| = HCP_LC P/All ratio for Jj 
C| = Cones per Jj tree 
Pj = Pollens per cone for Jj 

The number of trees of a Juniper species “i” per grid cell is calculated as 

Ti = GAPj x MODIS x TC 

Where 

GAPj = Fraction of Jj at 1 km grid (range 0-1) 

MODIS = MODIS derived percent tree cover per 1 km 2 grid cell (in fraction, range 0-1) 
TC = Tree count or number of trees. 


Ground truth (transect data) 


(a) Male to Female ratio 

(b) HCP_LCP to All ratio 

0 - Only enough cones to determine gender 

1 - Low Cone Production (LCP) tree 

2 - High Cone Production (HCP) tree 

(c) Cones per tree 

(d) Pollen grains per cone 
(d) Age 

(Height & edge effect) 


GAP derived distribution Pollen Source Mask (PREAM) 


Juniper Ashei 

December to January 
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Pollen Count at Source 
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GAP derived Juniper distribution 
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Juniper Pinchotii 

October-November 


Juniper 

Monosperma & Scopulorum 

March-April 
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GAP derived juniper distribution 
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Pollen Count at Source 
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California 
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GAP derived Juniper distribution 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 


Spatial resolution: ~1 km (990 m) 




PREAM 

Pollen Plume Simulation for Juniper Emissions 
For the period 15 December 2009 - 1 January 2010 

Run by Slobodan Nickovic, September 2012 


Atmosphere Model Setup) 
Model horizontal domain : Southwest US 
Model resolution: ~40 km 


Simulation period: 

15 December 2009 - 1 January 2010 


Boundary conditions: 1 degree global forecasts used to refresh 

• initial conditions every 24 hours 

• boundary conditions every 6 hours 



PREAM 

Pollen Plume Simulation for Juniper Emissions 
For the period 15 December 2009 - 1 January 2010 


PREAM_ - Pollen Regional/ \tmosjpherfc Model 
Derived from DREAM (dust), modified to simulate pollen 

- 4 particles bins 

- PREAM is online driven by the NCEP/ETA 
Emission: 

- Viscous-sublayer parameterization 

- Emission dependent on friction velocity 
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Model Validation ( Juniper pollen count/ m 2 ) 
Observed Forecast 


Sar>_Marcos,TX Dec-Jan 2009/2010 
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Report a Syndrome 


Demographics 


Syndrome Information 


Clinical Findings: Chronic Lung Disease Exacerbation 

-Symptoms (Reported by Patient) 




Productive Cough? 

O Yes 

ONo 

Nasal Discharge? 

OYes 

ONc 

Sore Throat? 

OYes 

No 

Wheezing? 

OYes 

O Nc 

Underlying Lung Disease (Asthma/COPD)? 

OYes 

; No 





Clinical Signs (from Physical Examination)- 


4 


Temp(C) O < 37.0 O 37.0 - 37.9 O 38.0 - 38.9 O 39.0 - 39.9 
Predominant Lung Findings 0 Rales O Wheezing Bilateral Unilateral 

Skin Rash? O Yes O No Oral Lesions? O Yes O No 
Lymphadenopathy? Yes No Diffuse Localized 

Splenomegaly? Q Yes O No Hepatomegaly? Yes O No 

Laboratory and X-Ray Data 




WBC Count: 0 < 5,000 O S,000 - 10,000 O 10,001 - 15,000 O > 15,000 

Platelet CL Q< 50,000 O 50,000 - 100,000 O 100.001 - 150,000 O > 15 
Chest X-Ray: 0 Normal O Abnormal 

M Infiltrate Q Hyperinflation l_l Cardiomegaiy Effusion 

GZ Sat. (Room Air) Q Normal (_ Abnormal 


Help Cancel Submit Report 
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National I Phenology Network 
Taking the Putse of Our Planet 

ABOUT PARTICIPATE RESOURCES 


EDUCATION 


RESULTS 
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Join Us! 


We are looking for volunteers to 
he] p us monitor plant and animal 
species found across the United 
Stales. Click "Observe" to Join us! 





ARCHIVE 


western columbine 

View All Species 


USA National Phenology Network 

The USA National Phenology Network brings together citizen scientists, government agencies, non-profit 
groups, educators and students of all ages to monitor the impacts of climate change on plants and animals in 
the United States, The network harnesses the power of people and the internet to collect and share 
information, providing researchers with far more data than they could collect alone. 


Learn more about us 


What is phenology? 


Phenology refers to recurring plant and animal life cycle stages, or phenophasestt, such as leafing and 
flowering, maturation of agricultural plants, emergence of insects, and migration of birds. Many of these events 
are sensitive to climatic variation and change, and are simple to observe and record. As an USA-NPN observer, 
you can help scientists identify and understand environmental trends so we can better adapt to climate change. 

Why is phenology important? 


USA-NPN News 


Phenology Feed Join the Conversation 


► Phen oclimatology Position at UA 

* Introducing the USA-NPN Video / 

^ Nature's Notebook: "How to Observe" Handbook ~ and Training Videos f 

► Phenology Special Issue in the Philosophical Transactions of the Royal Society 

* USA-NPN Reports (including Strategic Plan and 2009 Annua! Report)^ 

* Call for Papers: 4th Annual PROSE in Tucson, AZ, October 2010 ^ 


► Recem Media Reports 

► Newsletter Archive 


Sponsors 


^ Are you,,,? 

New to phenology? 
Ready to start observing? 
One of our partners? 
Interested in creating a partnership? 

An educator? 
Interested In finding data to use? 

A media outlet? 


About Ua | Contact Us \ Partners | Privacy | Aooeaaihifity I Glossary J Links \ Media Center 
USA National Phenology Network, National Coordinating Office 
IflSS E. Sixth SL, Tucson. AZGS721 
Admin Login 






signal to an extent that would allow detection, given accurate QA filtering and 
BRDF corrections. MODIS daily reflectances increased during the pollen season. 

S The DREAM model (PREAM) was sucessfully modified for use with pollen and 
may provide 24-36 hour running pollen forcasts. 

^Publicly available pollen forecasts are linked to general weather patterns and 
roughly-known species’ phenologies. These are too coarse for timely health 
interventions. PREAM addresses this key data gap so that targeting intervention 
measures can be determined temporally and geospatially. 

^The New Mexico Department of Health (NMDOH) as part of its Environmental 
Public Health Tracking Network (EPHTN) would use PREAM a tool for alerting the 
public in advance of pollen bursts to intervene and reduce the health impact on 
asthma populations at risk. 



